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An energy for the homogeneous cosmological models is presented. More specifically, using an 
appropriate natural prescription, we find the energy within any region with any gravitational source 
for a large class of gravity theories — namely those with a tetrad description — for all 9 Bianchi 
types. Our energy is given by the value of the Hamiltonian with homogeneous boundary conditions; 
this value vanishes for all regions in all Bianchi class A models, and it does not vanish for any 
class B model. This is so not only for Einstein's general relativity but, moreover, for the whole 
3-parameter class of tetrad-teleparallel theories. For the physically favored one parameter subclass, 
which includes the teleparallel equivalent of Einstein's theory as an important special case, the 
energy for all class B models is, contrary to expectation, negative. 
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I. INTRODUCTION 

Gravity is the only universal force; it is long range 
and dominates the cosmos. Energy has been one of the 
most useful physical concepts, no less so in gravitating 
systems — one need only recall its utility in the New- 
tonian Kepler problem. From the modern relativistic 
perspective energy-momentum is the source of gravity, 
yet — somewhat ironically — unlike all matter and other 
interaction fields in the absence of gravity, identifying a 
good description of the energy-momentum of gravitating 
systems has been the oldest and most controversial out- 
standing puzzle. Physically the fundamental difficulty 
can be understood as a consequence of Einstein's equiva- 
lence principle — from which it follows that gravity cannot 
be detected at a point. Hence the energy- momentum of 
gravity cannot have a proper local density; it is funda- 
mentally non-local and thus inherently non-tensorial and 
therefore non-covariant. (For a good discussion of this 
point see 

After Einstein proposed his gravitational energy- 
momentum density (see e.g. Q), various alternate pre- 
scriptions were introduced by other researchers (notably 
Papapetrou p|, Bergman and Thompson M0ller M, 
Goldberg @, Landau and Lifshitz Q, and Weinberg ||]). 
These investigations lead to a variety of expressions with 
no compelling criteria for favoring any particular one. 
Moreover these traditional energy-momentum pseudoten- 
sors are, as noted, necessarily not covariant objects, 
they inherently depend on the coordinates, so they can- 
not provide a truly physical local gravitational energy- 
momentum density. Caught between the equivalence 
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principle and the covariance principle, the pseudoten- 
sor approach has been largely questioned, although never 
completely abandoned. 

On the other hand, because the gravitational inter- 
action is local, some kind of local — or at least nearly 
local — description of the associated energy-momentum 
was still sought. The modern concept, introduced by 
Penrose to resolve this dilemma, is that properly 
energy-momentum is quasi-local, being associated with a 
closed surface bounding a region (for a nice review of the 
topic see 0). However it soon became clear that there is 
no unique quasi-local energy expression. Although there 
are some especially famous ones (e.g., Brown- York 
and Liu-Yau [IH) many definitions of quasi-local mass- 
energy have been proposed; they generally give distinct 
results. For example Bergqvist [13j studied several differ- 
ent quasi-local mass definitions for the Kerr and Reissner- 
Nordstrom spacetimes and came to the conclusion that 
not even two of the examined definitions gave the same 
result. 

Our view is that from the Hamiltonian perspective 
one can make sense of this situation — understanding not 
only why all these otherwise perplexing choices exist 
but also what is their real physical significance. Sim- 
ply put the energy of a gravitating system within a 
region — regarded as the value of the Hamiltonian for this 
system — naturally depends not only on the interior of the 
region but also on the boundary conditions imposed at 
the interface with the exterior (this is reasonable, after 
all the particular solution to the field equations depends 
on the boundary conditions). It has been found that 
the Hamiltonian necessarily includes an integral over the 
2-surface bounding the considered region. This Hamil- 
tonian boundary term plays two key roles: (i) it con- 
trols the value of the Hamiltonian, and (ii) its specific 
form (via the requirement that the boundary term in 
the variation in the Hamiltonian vanish) is directly re- 
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lated to the selected type and value of the boundary con- 
ditions. Many of the quasi-local proposals (namely all 
those which admit a Hamiltonian representation) can be 
understood in these terms: their differences are simply 
associated with different boundary conditions. Further- 
more, using a covariant Hamiltonian formalism, it has 
been shown [l5| that every energy-momentum pseu- 
dotensor can be associated with a particular Hamiltonian 
boundary term — which in turn determines the quasi-local 
energy-momentum that is respectively linked with the 
implied boundary conditions. In this sense, it has been 
said that the Hamiltonian quasi-local energy-momentum 
approach rehabilitates the pseudotensors and, moreover, 
dispels any doubts about the physical meaning of these 
energy-momentum complexes — since with this approach 
all their inherent ambiguities are given clear physical and 
geometric meanings. 

We want to emphasize that while there are many 
possible Hamiltonian boundary term energy-momentum 
quasi-local expressions — simply because there are many 
conceivable boundary conditions — it nevertheless has 
been found [l6j that in practice there usually is a par- 
ticular choice best suited to the task at hand. It has 
long been known that this is the case for familiar phys- 
ical systems where gravity is negligible. For example in 
thermodynamics we have not a unique energy but rather 
the internal energy, the enthalpy, and the Helmholtz and 
Gibbs free energies — each a real physical energy adapted 
to a specific interface between the physical system of in- 
terest and its surroundings; similarly in classical electro- 
statics the physically appropriate measure of energy ob- 
tained from the work-energy relation for a finite system 
depends on whether one considers fixed surface charge 
density or fixed potential — boundary condition choices 
which are respectively associated with the symmetric and 
the canonical energy-momentum tensor. 

Here we apply this Hamiltonian boundary term quasi- 
local energy-momentum approach to homogeneous cos- 
mologies. Our motivation is twofold. On the one hand 
the utility of having a good measure of the energy of such 
gravitating systems — a measure with sensible answers 
for the ideal exactly homogeneous case but which can 
be applied to any perturbations thereof and even quite 
generally — should be obvious. We are likewise motivated 
by the consideration that homogeneous cosmology affords 
an excellent set of models where one can test the suitabil- 
ity of our — or indeed any other — proposed (quasi-)local 
energy-momentum ideas. 

Accordingly it is important to here remark on our spe- 
cific choice of energy expression for these homogeneous 
cosmologies. We have noted five different approaches 
which lead us to exactly the same formula for energy- 
momentum for the homogeneous cosmologies in General 
Relativity (GR). (i) One may take the tetrad form of 
GR: among the many traditional energy-momentum ap- 
proaches M0ller's [l7| tetrad-teleparallel formulation has 
been highly regarded. This form of GR can be viewed 
from the gauge theory perspective, then a certain trans- 



lational gauge current (closely related to the energy- 
momentum expression proposed by M0ller) stands out 
(among other virtues it is the only classical expression 
that has an associated positive energy proof [18]). (ii) 
Another perspective is via the GR Hamiltonian: from 
the covariant Hamiltonian approach for GR one particu- 
lar Hamiltonian boundary term stands out as being the 
favored expression for general applications [l6| . In many 
situations, including the present application, it reduces to 
the tetrad-teleparallel gauge current expression, (iii) Al- 
ternately one can proceed from the spinor-parameterized 
Hamiltonian associated with the Witten positive energy 
proof [ltj [2(| [2l|. For an isolated asymptotically flat 
gravitating system the spinor field may be taken to satisfy 
Witten's equation, then the Hamiltonian has a positive 
value; however for a homogeneous cosmology the Wit- 
ten equation does not have appropriate solutions, instead 
it is appropriate to take the spinor field to be homoge- 
neous, which reduces the Hamiltonian to the aforemen- 
tioned expression, (iv) A fourth approach — which we like 
for its simplicity, generality and straightforwardness — 
is to consider the general teleparallel theory (which in- 
cludes Einstein's GR as a special case). For such the- 
ories, quite unlike the GR situation, investigators have 
advocated only two specific expression for the energy- 
momentum density (this is one of the virtues of treat- 
ing this whole class of theories rather than the one spe- 
cial case equivalent to GR), one of which — the aforemen- 
tioned tetrad-teleparallel gauge current — stands out as 
most suitable for this and most other applications, (v) In- 
stead one could begin directly with homogeneous cosmol- 
ogy: homogeneous cosmologies are naturally described in 
terms of a preferred homogeneous tetrad; then an energy- 
momentum expression based on the Hamiltonian for this 
preferred tetrad is clearly appropriate. 

There have been several studies aimed at finding the 
total energy of the expanding universe. An early inves- 
tigation proposed that our universe may have arisen as 
a quantum fluctuation of the vacuum. That model pre- 
dicted a universe which is homogeneous, isotropic and 
closed, and consists equally of matter and anti-matter. 
Albrow [221 ] and Tryon [23 ] proposed that our universe 
must have a zero net value for all conserved quantities 
and presented some arguments, using a Newtonian order 
of magnitude estimate, favoring the proposal that the 
net energy of our universe may indeed be zero. The gen- 
eral argument for this requirement is that energy can al- 
ways be represented by an integral over a closed 2-surface 
bounding the region of interest, so if the universe has an 
empty boundary then the energy should vanish (for a 
more detailed discussion of the vanishing of energy for a 
closed cosmology see [24[). Thus for closed cosmological 
models the total energy is necessarily zero. Years ago 
Misner [25| pointed out a technical problem with the at- 
tempts at explicit demonstrations of this statement. He 
noted that the integrands that were being used were ref- 
erence frame dependent (holonomic) pseudotensors and 
their associated superpotentials. None of those discus- 
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sions had specifically established exactly how these par- 
ticular non-tensorial objects behaved under changes of 
coordinates — but such changes were actually necessary 
in the calculations since the whole universe could not be 
covered with one coordinate patch. According to our un- 
derstanding it was Wallner [26| who first gave for GR a 
clear demonstration of this important vanishing energy 
for closed universe requirement; his integrand (effectively 
the same tetrad-teleparallel gauge current that we will 
use) was given in terms of a globally defined frame field. 

The subject of the total energy of expanding uni- 
verse models was re-opened by Cooperstock and Is- 
raeli (27l, Rosen [H|, Garecki [U, Jhori et al. [H, 
Feng and Duan (3lj , and others using various GR energy- 
momentum definitions. In one of these investigations 
the Einstein energy-momentum pseudotensor was used 
to represent the gravitational energy [28|, which led to 
the result that the total energy of a closed Friedman- 
Lemaitre-Robertson- Walker (FLRW) universe is zero. In 
another, the symmetric pseudotensor of Landau-Lifshitz 
was used [3(3] • Some works calculated the total energy of 
certain anisotropic Bianchi models using different pseu- 
dotensors, leading to similar results (32, [H, HI] . More 
recently Faraoni and Cooperstock [35| . in support of 
Cooperstock's proposal that gravitational energy van- 
ishes in the vacuum, argued (with the aid of a non- 
minimally coupled scalar source) that the open, or crit- 
ically open FLRW universes — as well as Bianchi models 
evolving into de Sitter spacetimes — also have zero total 
energy. Finally, calculations for the closed FLRW [36[ 
and some anisotropic Bianchi models [37| using the Ein- 
stein, Laudau-Lifshitz and other complexes in teleparallel 
gravity also led to vanishing energy. 

In the present work, we examine the energy- 
momentum (for both local and global regions) for a large 
class of cosmological models — more specifically we con- 
sider all 9 Bianchi types of homogeneous cosmological 
models using the Hamiltonian boundary term quasi-local 
approach in the context of not only GR but more gen- 
erally the tetrad-teleparallel theory of gravity 38]. This 
theory is a generalization of GR (it has been dubbed 
NGR, for new general relativity); a certain special case 
which is equivalent to Einstein's theory [lj, [3t| was 
first proposed by M0ller to solve the energy localiza- 
tion problem. This special case has been referred to 
as GR||, the teleparallel equivalent to Einstein's theory 
(a.k.a. TEGR); it has attracted the attention of several 
investigators, see e.g., pj, HI, HI 5H E3, El- Our mo- 
tive for this general approach is not only to improve our 
understanding of these cosmological models, but also to 
better understand the considered gravity theories and es- 
pecially to better understand the meaning and applica- 
tion of these quasi-local energy-momentum ideas. 

Most of the previous calculations of the energy of the 
universe used cartesian or spherical coordinate systems; 
here in contrast we use only the symmetry of spacetime 
given by the Bianchi group. In this major extension of 
previous results [36l l37| , by using the Hamiltonian for- 



malism with homogeneous boundary conditions we find 
that the energy — for all regions (any shape, large or 
small) — vanishes for all Bianchi class A models and does 
not vanish for any class B model. This result does not 
depend on the specific type of matter content of the uni- 
verse (e.g., baryons, radiation, cold dark matter, dark 
energy, cosmological constant); it depends only on the 
symmetry of the spacetime given by the structure con- 
stants of the Bianchi group. It is noteworthy that we 
find a simple physical difference related to the division 
of the Bianchi types into two classes. Moreover, for both 
GR and the physically viable teleparallel theories, our 
energy for Bianchi class B turns out to be, contrary to 
our expectations and presumptions, negative. 

The outline of the rest of this work is as follows: in 
the next section we discuss the Bianchi type universes. 
In Section 3 we briefly discuss the essential elements con- 
cerning the conserved energy-momentum expressions for 
Einstein's general relativity (GR), the general tetrad- 
teleparallel theory and the teleparallel equivalent of 
GR. In Section 4 we give the respective Hamiltonian- 
boundary-term quasi-local expressions. The energy cal- 
culation is presented in Section 5. This is followed by our 
concluding discussion of the remarkable results. 

II. BIANCHI TYPE UNIVERSES 

Suppose that the four-dimensional spacetime mani- 
fold can be foliated by homogeneous space-like hyper- 
surfaces of constant time t. Homogeneity means that 
each spatial hypersurface has a transitive group of isome- 
tries. The study of this 3-parameter isometry group, via 
the classification of 3-dimensional Lie algebras, led to 
the Bianchi classification of spatially homogeneous uni- 
verses [44j . 

These models are characterized by homogeneous — but 
generally anisotropic — spatial hypersurfaces parameter- 
ized by time. In a synchronous coordinate systems, in 
which the time axis is always normal to the hypersurfaces 
of homogeneity Et, we take the spacetime orthonormal 
(co) frame for these cosmological models to have the form 

?9°=dt, d a {t) = h a k {t)a k , a =1,2, 3 (1) 

where the three basis one-forms a k = a (x), k = 1,2,3, 
depend on spatial position in such a way that 

da k = X -C k lm a l A a m , (2) 

with C k i m being certain constants. The associated space- 
time metric then has the form 

ds 2 = -dt 2 +g lm {t)a\x)a m {x), (3) 

where gi m :— 6 a t,h a ih b m is a spatial 3-metric which de- 
pends only upon time; for our analysis it need not be 
diagonal. 
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There are 9 Bianchi types distinguished by the partic- 
ular form of the structure constants C k i m @, |45j|. They 
fall into two special classes: class A (types I, II, VIo, 
VII , VIII, IX) have A k := C m km = and class B (types 
III, IV, V, Vlft, Vll h ) are characterized by A k ^ 0. For 
our purposes here we hardly need any more details re- 
garding these types. We note that the respective scalar 
curvatures are vanishing for Type I, positive for Type 
IX, and negative for all the other types. Also, although 
the general idea of these models is homogeneous but non- 
isotropic, certain special cases can be isotropic; specifi- 
cally, isotropic Bianchi I, V, IX are, respectively, isomet- 
ric to the usual Friedmann-Lemaitre-Robertson- Walker 
(FLRW) models: k = 0, -1, +1. 

III. GR AND THE TETRAD-TELEPARALLEL 
THEORY 

Here we note some features, associated with conserved 
expressions for energy-momentum, for Einstein's general 
relativity (GR) and certain alternative tetrad-teleparallel 
theories, especially the one which is equivalent to GR. 

For all cases the analysis can be most conveniently ex- 
pressed in terms of differential forms and orthonormal 
(co)frames: 

d a = e a l dx\ i = 0,l,2,3. (4) 

In such frames the metric coefficients are constant: g^ v — 
diag(— 1, +1, +1, +1). We will consider on the same space 
more than one connection (equivalently, more than one 
covariant derivative or parallel transport rule). However, 
in all cases the connection is assumed to be metric com- 
patible: 

= dg^ - T x ^g\„ - T x u g^\ = -T^ - Y vli . (5) 

Hence the connection one-form coefficients are always an- 
tisymmetric: T a p = T\ a py Any particular connection is 
characterized by its curvature (2-form): 

R a p = dr a p + r% a rip, (6) 

and its torsion (2-form): 

T a = d$ a + T a pA^. (7) 

A. GR 

Einstein's general relativity is based on a Riemannian 
geometry, using the Levi-Civita connection — which has 
vanishing torsion. The Einstein field equations can be 
expressed in the form 

R a p A r] J '„ = -2kT m , (8) 

where ry Q/3 " := *(z9 Q A d 13 ■ ■ •) is the dual form basis. The 
expression on the left is the Einstein 3- form —2G L, IJ ,r] u , 



and the quantity on the right is the source energy- 
momentum 3- form T u ^rj U: with n = 87rG/c 4 . Here we 
take units such that c = 1. 

Let us first note that the above Einstein equation can 
naturally be rearranged into a certain special form (for 
similar arguments see 0, [47j ) : 

d^pAvJ^ = -T a pAdr } J tl -T a ^hT" l phr ] J lx -2KT tl . 

(9) 

The 2-form superpotential, T a p A rja 13 M , is the key ob- 
ject. Its differential is exact; consequently the 3-form on 
the right hand side is closed; hence it is automatically a 
conserved current — which includes the material energy- 
momentum linearly. The remaining pieces on the right 
hand side (which depend only on geometric quantities) 
can thus be interpreted as the energy-momentum density 
of the gravitational field. These forms arc just the quan- 
tities we wish to use to describe the energy of GR — and 
this rearrangement of the field equation is the easiest way 
to get them that we know of. 

Nevertheless, here we want to also present some ad- 
ditional material to (i) reinforce our thesis that these 
indeed are the most suitable energy-momentum expres- 
sions for homogeneous GR cosmologies and also to (ii) 
include a large class of alternate theories. 

Hilbert noted that the above Einstein equation |5J can 
be obtained from the scalar curvature Lagrangian den- 
sity by regarding it as a function of the metric. Later 
certain variations on the theme were considered. For 
the vacuum — or even non-derivative coupled sources like 
scalar fields or the Maxwell and Yang-Mills gauge fields — 
one can simply vary the coframe and the metric compat- 
ible connection one-form independently Then the con- 
nection variation yields an equation which requires the 
connection to have vanishing torsion, i.e., to be the Levi- 
Civita connection of GR. Let us next note a result [48[ we 
can get most simply by proceeding from that Lagrangian 
formulation. 

The scalar curvature Lagrangian 4-form, Cb, = Rf] = 
R af3 A r] a p, can be rearranged as indicated: 

R°p A r)J = (dT a p + T% A P» A nj = 

d (r a p a V J) + p> a d v J + r% a n> a r,J. (10) 

Now by dropping the total differential we obtain 

2 K £ tctl = -dzT A T a p A r,J^ + T Q 7 A V p A r)J, (11) 

a modified Lagrangian density which gives the same field 
equations. It is noteworthy that the canonical momen- 
tum conjugate to the co-frame, 

d£tcti _ ra p _ TOl p , TP7 1 

is just the already encountered superpotential 2-form. 

Recall that the Levi-Civita connection is linear in the 
differential of the frame: 

= \ [(d# a )fr + (d#p) ja + [d^)p a ] W . (13) 
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This relation can be used to eliminate T a p from (fTTj) to 
obtain a certain specific action: Z^GRtet = ^GRtet ($> d$) , 
which is quadratic in di9. In addition to this tetrad ver- 
sion of Einstein's theory there are other more general 
tetrad theories of gravity which follow from other La- 
grangians quadratic in Ad. Such tetrad theories are some- 
what interesting in their own right. Let us consider them 
in the next subsection. 



B. tetrad-teleparallel theory 

In this work we are concerned especially with the 
general homogeneous Bianchi cosmologies. Such mod- 
els have a natural preferred global homogeous orthonor- 
mal frame (a.k.a. tetrad, vierbein). If a geometry has a 
preferred tetrad one can naturally introduce a new par- 
allel transport rule (i.e., a new connection) such that 
this frame is, by definition, parallel. A geometry with 
a preferred global parallelism is referred to as teleparallel 
(a.k.a. absolute parallel). Hence, especially for the homo- 
geneous cosmologies, an appropriate theoretical geomet- 
rical framework is the tetrad/ teleparallel formulation. 

Conversely, for any teleparallel geometry by definition 
the curvature vanishes and the parallel transport is path 
independent. It follows that a global preferred frame may 
be constructed by starting with any orthonormal frame 
at a single point and parallel transporting it along any 
path to all the other points. The resultant global tetrad 
field is unique (up to an overall global constant Lorentz 
transformation). In this constructed frame the connec- 
tion coefficients vanish. Hence in a teleparallel geometry 
there is (up to a an overall global constant rotation) a 
preferred frame in which the connection also vanishes. 
The basic variable can then just be taken to be this pre- 
ferred tetrad, which is most conveniently represented as 
the (co)frame i? Q = e a idx % . Since the connection co- 
efficients vanish in this frame the associated teleparallel 
torsion 2-form is then simply given by the frame differ- 
ential: 



T^ al3 d a A d p = M a 



(14) 



Thus wc can forgo further mention of the teleparallel con- 
nection here and simply regard a teleparallel theory as a 
theory for a preferred tetrad. 

Tetrad theory field equations can be obtained from a 
Lagrangian 4-form £ tot = £ tct {tf a ,d$ a ) + £ mat (i?). The 
variation with respect to the frame gives 



SCt 



W A 



ddW 



•<50" A 



fdCt 



(15) 



This expression has the form 

<5Aot = Sd^ A r M + <5iF A (-t M - %) = 

d (60* A t„) + 6&* A (dr M - t„ - TJ, (16) 



which identifies as the canonical conjugate field mo- 
mentum and i M and 7^ as, respectively, the 3-forms of 



gravitational energy-momentum and the material source 
energy-momentum density. Using (|16p and Hamilton's 
principle gives the tetrad field equation in a form rem- 
iniscent of that used by Einstein in his search for his 
gravity field equations [49( : 



dr., = t L 



(17) 



The right hand side is naturally a conserved current 3- 
form: 



d(t M + = 0, 



(18) 



(i.e., because the lhs is exact the rhs must be closed). 
The associated conserved total energy-momentum within 
a volume V is 



Pm(V)= t ll + Tp= d V = f T(l . (19) 

JV JV JdV 

The tetrad-teleparallel formulation is natural from the 
gauge theory point of view. From that perspective the 
above expressions are those of the translational gauge 
current (the conserved current associated with spacetime 
translations according to Noether's first theorem, see, 

e.g., EE Gil 53, EE El E3). 

For the tetrad-teleparallel theories the canonically con- 
jugate momentum field, 

dC 1 



a/3 



(20) 



ii ■ dTtl 2 i* 

is generally taken (in order to have quasi-linear second 
order field equations) to be a linear combination, 

t = «T (<2iT ten + a 2 T vcc + a 3 r axi ), (21) 

of the tensor, vector, and axivector irreducible parts of 
the teleparallel torsion [50] : 



Ti a 
vcc fj,u 



T' 



/z]A) 



rpCtflV 

axi 



(22) 

There is thus a 3-parameter class of such theories [38| . 
The generic theory determines a preferred frame. It 
turns out, however, that one special parameter choice 
(4a3 = a2 = — 2ai) is distinguished: it actually has lo- 
cal Lorentz gauge freedom. This model, with a\ = — 1, 
has been encountered hereinbefore (fT2l): it is known as 



GRtet or GR||, the teleparallel equivalent of Einstein's 
GR (a.k.a. TEGR), and was first proposed by M0ller [13 
to solve the GR energy localization problem. 



IV. QUASI-LOCAL BOUNDARY EXPRESSION 

As mentioned in the introduction, for Einstein's 
GR many energy-momentum expressions — both quasi- 
local and reference frame dependent pseudotensors — 
have been proposed. It should be emphasized that, de- 
spite much effort and many nice results, there is no con- 
sensus as to which, if any, is best. The Hamiltonian ap- 
proach certainly helps. From that perspective the energy- 
momentum is determined by the boundary term in the 
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Hamiltonian [Tj, [5l|]. Although (at least formally) the 
formalism allows for an infinite number of Hamiltonian 
boundary expressions (including all the superpotentials 
that generate the pseudotensors), the ambiguities have 
been tamed: each expression has a geometrically and 
physically clear significance associated with the boundary 
conditions determined from the variation of the Hamil- 
tonian 47] . Nevertheless there are (at least formally) an 
infinite number of possible boundary conditions. 



A. Tetrad-teleparallel energy-momentum 

On the other hand, the situation for the 
tetrad/teleparallel theory is in sharp contrast. In- 
vestigators [HI p, M, m, M, S3, S3, S3, M, HI were 

led to only two [5J] closely related quasi-local boundary 
term expressions for the energy-momentum within a 
volume V: 



Pi{V) 



Jav Jov 



(23) 



they are, respectively, the Ivfoller 1961 [TtJ expression 
and the translational gauge current derived above. M0ller 
had pointed out that his superpotential (which appears 
here as a 2-form integrand) is tensorial (i.e., it transforms 
homogeneously under a change of coordinates) ; however 
its differential, 



d(e" 



de^, A t„ 



; dr., 



(24) 



the M0ller tetrad-teleparallel energy-momentum 3-form, 
is not a tensor with respect to coordinate transformations 
(as M0ller himself noted) — because of the factor de'V In 
contrast, it should be emphasized that both the transla- 
tion gauge current superpotential 2-form and its dif- 
ferential, the gauge current 3-form (|17[) . are true tensors 
under changes of coordinates. Generically, the tetrad- 
teleparallel theory has a natural preferred frame (no local 
frame gauge freedom), then the translational gauge cur- 
rent energy-momentum expressions have no ambiguity at 
all. 

However for the one special case of greatest interest, 
GRtet, the theory does have local Lorentz gauge free- 
dom. In that case the gauge current expressions do de- 
pend on the choice of orthonormal frame, and thus still 
contain some observer dependent information mixed in 
with the physical information in the energy-momentum 
expression. Nevertheless we can regard the gauge cur- 
rent expressions as preferable to any of the pseudoten- 
sors or Moller's 1961 expression, since dependence on an 
orthonormal frame is more physical than dependance on 
an arbitrary choice of coordinates. 

Concerning the ambiguity of the choice of frame for 
GRtet, it is important to note that the quasi-local val- 
ues depend only on the choice of frame on the boundary, 
and not on the choice within the interior of the region. 



Moreover in the case of interest here (homogeneous cos- 
mologies) there is a preferred frame — and thus there is 
no ambiguity at all. 

In summary, in the tetrad-teleparallel formulation un- 
like GR there is no big ambiguity in the choice of 
expression. Thus our consideration of this general 
tetrad/teleparallel class of theories yields two benefits: 
(i) it allows us to get a result of great generality, apply- 
ing to this whole class of theories, (ii) moreover, when 
specialized to TEGR, it determines a unique preferred 
energy-momentum expression for GR. This expression, 
the tetrad-teleparallel gauge current has been regarded 
as one of the best, perhaps the best, description of the 
gravitational energy-momentum for GR. 



B. The covariant Hamiltonian approach 

On the other hand, let us also briefly discuss the al- 
ternative of taking the usual Ricmannian geometry ap- 
proach to GR. A covariant Hamiltonian formulation has 
been developed for general geometric gravity theories. 
The Hamiltonian which generates the evolution of a spa- 
tial region V along the vector TV is given by an integral 
of the form 



H{N, V) = I N^H 



ov 



B(N). 



(25) 



It turns out that H,^ is proportional to field equations 
and thus has vanishing value (e.g., for GR ([5]) we have 
= ~{2n)~ 1 R a i3 A nji^ - %). The Hamiltonian 
boundary term controls the Hamiltonian value and the 
boundary conditions. It has considerable freedom. The 
analysis led to certain specific quasi-local Hamiltonian 
boundary term 2-form expressions related to various 
types of boundary condition choices (e.g., Dirichlet, Neu- 
mann) [l5l . [47], [5lT | for quite general gravity theories. 
When specialized to Riemannian GR, it was noted that 
one of these expressions is singled out, in that it corre- 
sponds to boundary conditions imposed on a complete 4- 
covariant object (the tetrad) and gives the desired Bondi 
energy flux [JJ|: 



2kB{N) = AT a A i N nJ + D N a Ar] a 



(3 



(26) 



Here N describes a spacetime vector field which selects 
the components of energy and momentum, Ar := T — 
r and An := n — fj, where T and 77 are reference (or 
ground state) values. (We note that the same energy- 
momentum expression in its holonomic form was found 
by Katz, Bicak, and Lynden-Bell. They have extolled its 
virtues in several works, see in particular (55| and the 
references therein.) 

Appropriate choices for the Homogeneous cosmologies 
are to take the frame to be the preferred Bianchi frame, 
the reference connection to have vanishing components 
in this frame, and for the spacetime displacement vector 
field N to be homogeneous, i.e., to have constant com- 
ponents in this frame. Consequently the second term 
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vanishes and the boundary expression reduces to 

2kB(N) = T^A i NVa0 = N^^lJ^rp. (27) 

This is, as promised, the same superpotential form en- 
countered earlier (|9ll2p . an expression whose utility has 
been recognized at least since [191 ] . This succinct argu- 
ment shows that the preferred Riemannian GR covariant 
Hamiltonian quasi-local boundary term coincides in this 
situation with the GR|| gauge current. The same expres- 
sion also follows naturally from the Hamiltonian formula- 
tion of the tetrad-teleparallel theory [l8| and, moreover, 
from the Hamiltonian boundar y te rm associated with the 
Witten positive energy proof fill I20I l2l| with homoge- 
neous spinor field. 

V. THE BIANCHI ENERGY CALCULATION 

As we have argued, one can regard GR as a special case 
of the tetrad-teleparallel theory, so we consider in detail 
the latter formulation. The energy-momentum integral 
over the boundary of a region at a fixed time t is 

P»(V) ■= I V = / (28) 

JdV JdV 1 

In this Bianchi cosmology case the components of t /j q/3 
and T^ a p (in the preferred teleparallcl frame) are func- 
tions of time alone — dependence on the spatial coordi- 
nates shows up only in the teleparallel coframe via 
the a m . Consequently, in detail 

P„(V) = \r^{t) I Va p = \r, a0 (t) [ d Vafi (29) 

z JdV z JV 

= ^(t)[T% + T\ a S$ - T\ S° a ](t)V. 

Note that the energy-momentum is given by an integral 
over the 2-dimensional boundary of the region, yet it 
turns out to be simply proportional to the size of the 
included 3-dimensional volume (not its shape); it is also 
noteworthy that there is no dependence on the location 
(this is just as it should be considering homogeneity). 
Effectively, because of the homogeneity, we get a homo- 
geneous energy-momentum density, and a unique local- 
ization of energy-momentum. 

To properly appreciate our results here, one should 
note that most previous calculations of energy for cos- 
mological models considered the isotropic FLRW mod- 
els, or one or at most a few Bianchi types (and they 
then usually confined their results to diagonal metrics). 
They generally used holonomic energy-momentum ex- 
pressions, which correspond to rather different boundary 



conditions. Typically they considered some ball of con- 
stant radius around the origin or only the whole space. 
The expressions were not manifestly homogenous; the re- 
sults were shape dependent, not simply proportional to 
volume. 

Continuing with our calculation, we need to find the 
explicit value of the teleparallcl torsion (fl"4"l) . From ((T|) 
we find 

T° = 0, T a = d{h a k a k ) = h a k dtAa k + h a k da k . (30) 

using this along with @ we obtain the torsion tensor 
components T ^ = and 

T a ob = h a k h k bl T a bc = h a k C k i m h l b h m c . (31) 

From the latter it follows that 

T 7 & 7 = T c bc = C m kmh k b = A k h k b . (32) 

Consequently the energy-momentum is 

P,(V) = T^{t)T^,{t)8lV (33) 
= T, ob (t)T c bc (t)V = T^ ob (t)A k h k b (t)V, 

which vanishes for all class A models. (The physical in- 
terpretation is that the negative gravitational binding en- 
ergy density exactly cancels the positive material energy 
density.) 

For class B models we need 

KTfj° b ~ (aiTtcn + G^Tvcc + 03?axi) fj,° 

= [oiT + (o 2 - ai)T vcc + (a 3 - ai)T axi ] M ob .(34) 

For energy we need just the fj, = component. As T ^ — 
and T ax ; is totally antisymmetric it reduces to 

nr m = ( a2 -ai)T voc o 0fc = (a 2 - ai )iT c c h 

= (a 2 - ai )^A k h k a S ab . (35) 
Hence the energy within any volume V at time t is 

P°(V) = g m r 0b (t)A k h k b (t)V = ^-^A k A l9 k \t)V. 

6k 

(36) 

From various investigations it has been found that cer- 
tain parameter restrictions should be imposed. In par- 
ticular for the proper Newtonian limit we must require 
the so-called viable condition [50| . 2ai +02 =0. More- 
over normalization to the Newtonian limit gives a\ = —1. 
This leads us to the one parameter teleparallel theory, 
also known as NGR (New General Relativity (38|). If 
these parameter conditions are satisfied we have, for all 
viable cases (including the very special case of GR|| , for 
which we also have a\ + 2aj, = 0), 

E = P°(V) = -K- l A k A l9 k \t)V{t) < 0, (37) 
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i.e., negative energy for all regions in class B mod- 
els. With similar computations these same energy re- 
sults can be directly verified for GR using the relations 
(I12I13I19I27I) . 

Similarly one could find the explicit value of the "linear 
momentum" P C (V) [5y|. However, to get this one needs 
to evaluate 

T c ob = k- 1 [aiT + (a 2 - ai)T vcc + (a 3 - ai)T axi ] c ob , 

(38) 

which is linear (but not so simple) in T c o& = h c jh? V The 
calculation is straightforward but not enlightening. 



VI. DISCUSSION AND CONCLUSIONS 

We have obtained some new insight regarding the en- 
ergy of homogeneous cosmologies, especially from con- 
sidering the questions: what energy should be associated 
with a region of the universe? Does the total energy of a 
closed universe vanish? Is the energy of an open universe 
positive? Does the energy of empty flat space vanish? 

Specifically, using a natural prescription we found the 
value of the (quasi-)local energy-momentum for the gen- 
eral tetrad-teleparallel theory (which includes Einstein's 
GR as an important special case) for all 9 Bianchi types — 
with general homogeneous gravitational sources. From 
the Hamiltonian approach we found that, for comov- 
ing observers with homogeneous boundary conditions and 
reference, the energy vanishes for all regions in all Bianchi 
class A models, and it does not vanish for any class B 
model. This is the case for the whole 3-parameter class 
of tetrad-teleparallel theories. According to our measure, 
the one parameter set of viable teleparallel theories with 
a good Newtonian limit, which includes the teleparallel 
equivalent of Einstein's theory, has negative energy for 
all class B models. 

We note that all the cosmologies in the Bianchi class A 
models can be compactified, so our vanishing energy re- 
sult is consistent with the requirement of vanishing total 
energy for a closed universe. The class B models can- 
not in general be compactified, however there are some 
special exceptions. Nevertheless these exceptions are not 
counter-examples: while in certain cases the metric ge- 
ometry can be compactified, this cannot be done in such 
a way that the frames match up to give a globally defined 
smooth frame (i.e., they are not globally homogeneous). 
Indeed a proof of this used a calculation virtually 
identical to ours above without noting the energy inter- 
pretation. 

It is noteworthy that our energy depends only on the 
symmetry of the spacetime given by the structure con- 
stants of the Bianchi group, and that our result holds for 
all types of material sources including dark matter and 
dark energy (either as a cosmological constant or as some 
kind of unusual field like quintessence appearing as a part 
of the energy- momentum tensor). 

The results presented here can be specialized to the 



few Bianchi models which can be isotropic. For those 
cases, our homogenous results have been compared with 
those found using a similar approach for the more familiar 
"isotropic-about-one-point" FLRW formulations [H, H^] . 
The isotropic Bianchi I is isometric to the flat k = 
model, and both are found to have vanishing energy. The 
isotropic Bianchi IX is isometric to the k = +1 model; 
the energy of the latter for spherically symmetric regions 
first increases and then decreases to zero as the radius 
reaches the antipode so that the volume encompasses the 
whole close space (thereby satisfying the vanishing en- 
ergy for any closed space criteria), whereas the Bianchi 
IX has zero energy for all regions. The open k = — 1 
FLRW model was found to have negative energy, quali- 
tatively but not quantitatively like its isometric Bianchi 
counterparts V and VII^. It is remarkable that one spe- 
cial case, with scale factor a(t) = t, is actually isometric 
to empty Minkowski space with the expanding spatial 
slicing, t — \fT 2 — R 2 . The energy is negative for our 
measure applied to the expanding space but, of course, 
the energy vanishes for Minkowski with the constant T 
slicing. Thus the open and closed models offer differ- 
ent "localizations of energy" for the same physical situa- 
tion, and provide good examples of the effect of different 
boundary conditions and time evolution vectors. 

Although our energy values are obtained as integrals 
over the boundary of a region, they are not truly quasi- 
local. Properly a quasi-local quantity depends only on the 
physical data on the boundary of a region. Our values 
are obtained from integrals of the Hamiltonian boundary 
term over the boundary, but the boundary values, the ref- 
erence values, and the displacement vector N are selected 
using the global homogeneity. This gives us homogeneous 
localizations. For these models the homogeneity gives a 
physically meaningful preferred localization. 

Suppose one were to take a region of a homogeneous 
cosmology and regard it as an isolated gravitating sys- 
tem with the exterior being empty, then far away the 
isolated system would act like a Newtonian mass point 
and the asymptotic gravitating field should be nearly 
static. Then there are strong fundamental arguments 
(see [111 [Ml [sH, |H, HH) that the total effective mass- 
energy — physically determined by the parameters of a 
large Kepler orbit but mathematically defined as a cer- 
tain 2-surface integral at spatial infinity — should be pos- 
itive (i.e., gravity is attractive), and it should van- 
ish iff the system is empty Minkowski space. In par- 
ticular an isolated gravitating system would allow en- 
ergy to be extracted — or if left alone it would sponta- 
neously radiate — until it reached the lowest available en- 
ergy state. By scaling the lowest energy cannot have a 
finite negative value. A lower bound of negative infin- 
ity allows one to extract an infinite amount of energy, a 
source of perpetual motion, contrary to a fundamental 
thermodynamics principle — such a system would never 
be stable, it would continue to radiate forever. Thus the 
lowest state should be non-negative. The ground state of 
geometric gravity is Minkowski space with zero energy. 
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From this perspective it is reasonable to expect that 
the quasi-local energy determined by a suitable boundary 
integral at a finite radius, would also be non-negative and 
would vanish iff the interior were flat Minkowski space; 
indeed these properties have been regarded as desiderata 
for any good quasi-local energy (see, e.g., [IoLl65l|). 

The measure of energy used here for homogeneous cos- 
mologies is based on a different concept. Here we regard 
our cosmological region not as an isolated system but 
rather as existing in an exterior which is just like the 
interior. This is the essential significance of our choice 
of homogeneous boundary conditions, reference and evo- 
lution. One consequence is that our cosmological en- 
ergies do not satisfy the two important aforementioned 
quasi-local desiderata: indeed for the expressions con- 
sidered positivity need not hold, and zero energy iff flat 
Minkowski space does not hold — in either direction. 

Regarding E — 0. For any homogeneous measure of en- 
ergy in Bianchi class A models it is certainly quite reason- 
able to have a vanishing value for all regions, since a ho- 
mogeneous energy density energy must necessarily vanish 
at least for all compactifiable regions — and these models 
can be compactified (this is most obvious for Bianchi I 
with 3 torus topology identifications): the energy of a 
closed universe must vanish. This does not conform to 
the standard quasi-local criterion of a unique quasi-local 
E = Minkowski ground state. 

The most remarkable feature of our results is that 
many cosmological models were discovered to have neg- 
ative energy. It had been expected that we would al- 
ways find non-negative energy. One consequence of this 
strong belief was that a minor sign error was, quite un- 
fortunately, overlooked in some calculations; this lead to 
an incorrect claim of positive energy for Bianchi B mod- 
els bein g re ported in some recent conference proceed- 
ings [581 l59j|. 

It should be noted that a negative energy value for 
these open Bianchi cosmologies is not just a peculiar fea- 
ture due to our choice of expression or model. Indeed the 
results of calculations (reported in the aforementioned 
conference proceedings), done using essentially the same 
principles articulated here, likewise led to negative energy 
for the open FLRW (homogenous and isotropic) model. 
Furthermore, according to our simple calculations, most 
commonly used energy expressions will give the same en- 
ergy signature for both the open FLRW and Bianchi 
cosmologies. Yet, remarkably, we have not found any 
conspicuous reports of this simple fact, although there 
is some evidence which implicitly supports our energy 
signature. In particular Bannerjee and Sen [32| evalu- 
ated the Einstein pseudotensor for some Bianchi mod- 
els. Among other results, they reported a non-vanishing 
energy for finite regions in some class B models (specifi- 
cally, types III, V, VI), but they did not remark that the 
sign (easily evaluated) of their calculated energy value for 
these models is actually negative. The Landau-Lifshitz 
pseudotensor was used by Jhori et al. [3(| to calculate 
the energy for the closed k = +1 FLRW model; one can 



see from their Eq. (10) that they would have found a 
negative value for the k = — 1 model. 

The stability argument for non-negative energy is 
quite compelling — for isolated gravitating systems, which 
should be settling down to a stable equilibrium state. 
The cosmological models considered here, however, are 
quite different in kind. They are inherently dynamic, 
very unlike stable isolated systems. We do not see that 
there is any fundamental objection to assigning to them 
a negative energy value. 

Moreover, according to the following argument it is 
quite reasonable that negative spatial curvature models 
may have negative energy. Consider that a region of pos- 
itive spatial curvature is like a convex gravitational lens, 
tending to focus light rays, acting as if it had a posi- 
tive matter density which attracted the rays. Whereas a 
negative spatial curvature region acts like a concave lens, 
causing light rays to be defocused, acting as if it had a 
negative matter density which repelled light rays. Our 
definition of energy is actually the value of the Hamilto- 
nian but this can be expected to have some correlation 
with the effective active gravitational mass of the region. 

Identifying a good measure for the energy of a grav- 
itating system has remained an outstanding problem. 
Here we considered this issue for any region in homo- 
geneous cosmologies. We presented for all models the 
energy value associated with appropriate boundary con- 
ditions as given by our favored covariant Hamiltonian- 
boundary-term approach. We discovered, in contrast to 
the expectations, zero energy for some non-flat cases and, 
most surprisingly, negative energy for many (all Bianchi 
Class B). 

Finally, let us again emphasize that our specific re- 
sult hinges on the particular chosen measure of energy- 
momentum. We noted five perspectives which led to 
the same energy-momentum expression for GR: certainly 
the chosen expression is a natural one for both homoge- 
nous cosmologies and the teleparallel theories, and it co- 
incides with the GR|| tetrad-teleparallel gauge current. 
Yet our most fundamental argument for distinguishing 
this energy-momentum expression is via the covariant 
Hamiltonian approach — which includes and goes beyond 
the Noether analysis [13] • It should be noted that we ap- 
ply the Hamiltonian analysis to the general theory, and 
then impose the homogeneous symmetry on the result- 
ing expressions. It has long been known that one cannot 
in general first impose the symmetry and then do the 
Hamiltonian analysis. That approach happens to work 
successfully for Bianchi class A models but does not work 
for class B models [57]. 

Our energy-momentum is just the value of the Hamil- 
tonian (which dynamically generates the space-time) 
with homogeneous boundary conditions, reference and 
time slicing. We emphasize that our approach is entirely 
geometric (i.e., coordinate independent). Certainly there 
are other measures of energy-momentum — and they may 
have their own particular virtues — but we do not see how 
there can be any other which is more suitable for all 
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the homogeneous cosmologies. The results reported here 
can be considered as a standard measure of the energy- 
momentum for homogeneous cosmologies. 
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